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Abstract 


Radiative taps can be made in the cladding of optical fibers using the process of laser 
etching. Study is being carried out on the effect of such taps on the dispersion charcteristics on 
an equivalent planar dielectric waveguide. A hybrid spectral/mode-matching method has been 
adopted to calculate the electromagnetic fields guided by such a structure. The advantage 
of the spectral method is that the Fourier transform automatically enforces the boundary 
conditions on the taps, and the mode-matching method is straight forward and the solution 
lends itself to physical interpretation easily. The method thus combines advantages of both 
the methods. The results show a strong dependance of tap parameters on the dispersion 
characteristics. This can be utilised in the making of dispersion shifted and dispersion flattened 
fibers. 
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Chapter 1 
Introduction 


Traditional means of carrying out signal processing in optical communication 
is to convert the optical signals into electrical signals using photo detectors, 
and then carry out the processing electronically. However in such a scheme 
the speed of operation is limited by the speed of operation of the associated 
electronic circuits, which results in under utilisation of the enormous band- 
width that a optical fiber can offer. Thus, various innovative methods are 
being put into use to carry out the processing of signals at optical frequen- 
cies. A majority of the methods adopt one or the other means to tap out 
signals from a fiber without causing any major disturbance to the modes 
propagating inside the fiber. If a series of taps are distributed on an optical 
fiber, we obtain a fiber optic tapped delayline. 

Figure 1.1 shows a typical fiber optic delayline filter. Optical energy is 
radiated out of the taps made in the cladding of the fiber. Each beam is 


] Detector 



Figure 1.1: Typical fiber optic delayline filter 


made to pass through a light modulator, each beam thus passes through 
materials of different transmittivity. The beams are then collected by a lens 
and focussed on a detector. Depending on the transmittivities and distance 
between taps, the detector output changes and thus, the delayline acts as a 
signal processing filter. 

1.1 Methods to tap out optical signals from 
fiber 

Several approaches [2] - [4] , [13] to fiber tapping have been proposed such 
as macro-bending, polishing of the cladding for evanescent couping, slicing 
and deposition of dielectric mirrors in fiber core, and epoxy attachment of 
another fiber, local heating to induce divitrification and hence scattering, 
etc. All the methods mentioned above involve mechanical techniques and 
thus are not easy, inexpensive and reliable. 

A new technique for fabrication of high quality micron-scale fiber-optic 
taps using ablative chemical etching (LACE) was introduced by Imen et al 
[2], This technique involves the use of a CO 2 laser to etch a cut into the 
cladding of the fiber. The fields propagating in the cladding are refracted 
and reflected by the facets of the cut thus giving us radiating beams while 
not significantly affecting the signal propagating in the core. 

The set up for the fabrication of taps on optical fibers using the LACE 
is shown in fig 1.2 . A high power CO 2 laser beam is focussed onto the fiber 
using a lens and a tap is made by repeated translation of the beam across 
the fiber. The whole set-up is kept in a reactive gas ambient so as to stop the 
formation of debris at the edges. As the laser beam has a gaussian profile 
the resultant cut is gaussian in shape. By controlling the power, pulse time, 
polarisation and spot size of the beam, the scanning speed and the number 
of scans a high degree of control is achieved over the depth and width of the 
cut. With the lateral movement of the fiber along with laser scanning we can 
obtain shapes other than gaussian, viz ., triangular, rectangular etc. 

Radiation properties of such taps have been studied using the traditional 




Figure 1.2: The LACE fabrication process 


approach representing the cladding fields as rays at various angles by Imen 
et al [2] . However, it is also important to study the effect of a distribution of 
LACE taps on the dispersion characteristics of a fiber. A complete analysis 
of the waveguide dispersion characteristics of a fiber with LACE taps, for 
various tap geometries may open up applications of tapped fibers in the 
making of dispersion flattened fibers, dispersion shifted fibers etc. 

1.2 Problem overview 

In an optical communication system the bandwidth is limited mainly by the 
spectral width of the source and the dispersion characteristics of the opti- 
cal fibers. Dispersion in optical fibers can be of various types [17]. There 
is material dispersion, which occurs because of the fact that the refractive 
index varies as a nonlinear function of the optical wavelength. As a result 
various spectral components of a given mode travel at different speeds, de- 
pending on the wavelength. Waveguide dispersion, which is dependent on 
the fiber configuration. Depending on the waveguide configuration the prop- 
agation factor is a nonlinear function of the frequency and hence results in 
different speeds for different optical frequencies. Another form of dispersion 
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Figure 1.3: Waveguide dispersion characteristics of optical fiber 


is the intermodal dispersion, which is a result of different values of group 
delay for each mode at the same frequency. This mechanism does not af- 
fect single-mode operation. In this thesis we investigate the effect of finite 
number of LACE taps distributed periodically on the waveguide dispersion 
characteristics of a single mode fiber. 

In this thesis analysis has been carried out for the case of a rectangular 
slab waveguide with taps instead of optical fiber, since it has been shown that 
the intensity of fundamental TE mode of a rectangular slab is almost similar 
to LPqi mode of a weakly guiding optical fiber for the same core width and 
refractive indices of core and cladding [17], [14], [10]. For the case of a weakly 
guiding fiber the difference in the core and cladding refractive indices is very 
small. In such a fiber the expressions for the field components become simple 
and the fields thus obtained are similar to the fields in a rectangular slab with 
the same refractive indices, and with the width of core taken to be equal to 
the diameter of the core of the fiber. 

A hybrid spectral/mode- matching method has been adopted to solve the 
thesis problem [1],[5]. The method does not employ any approximations 
and can be applied practically to any dielectric waveguide. Moreover, the 
solution lends itself to physical interpretations very easily. The method also 
• J -- -i — *■ mutual couDline between different modes, which 



are generated due to the presence of taps. Other methods like FDTD [15], [9] 
were looked at, but the FDTD method requires at least 12 samples per 
wavelength which in optical domain gives around 0.1 micron, and also from 
the fact that the total fields are computed in every iteration. That would 
result in prohibitively large computation time when dealing with the nature 
of taps and, distances we were looking at. 

The method we use, calculates the electromagnetic fields in different re- 
gions of the structure. The solution is then constructed by choosing a suitable 
combination of these waves so as to satisfy the boundary conditions. The 
advantage of the spectral method is that the Fourier transform automatically 
enforces the boundary conditions on the taps. The method thus results in 
efficient calculation of the electromagnetic fields guided by the structure. 

1.3 Overview of the results obtained 

From the results of the simulation that have been carried out it has been 
observed that the propagation constant increases when taps are made in the 
cladding of the fiber. The propagation constant increases as the number 
of taps was increased. The propagation constant increases with an increase 
in tap depth. The increase in tap width does not cause much change in 
propagation constant. The results show that the propagation constant at a 
particular wavelength can be changed considerably by varying various tap 
parameters. 

1.4 Organisation of the thesis 

In chapter 2 we present the detailed discussion of the spectral/mode matching 
method. In chapter 3 we discuss the simulation package that has been devel- 
oped based on the spectral/mode matching method. In chapter 4 we present 
the results of simulation experiments that have been carried out for various 
tap parameters using the software that has been developed. In chapter 5 we 
present the conclusions and scope for further work. 



Chapter 2 


Spectral /Mode-matching 
method 


In this chapter we discuss the hybrid spectral/mode- matching technique, that 
has been adopted to solve our thesis problem. The method though rigorous 
in approach results in efficient calculation of the electromagnetic fields guided 
by the tapped optical fibers. The strategy adopted is to calculate fields in 
different regions of the structure. The boundary conditions are then applied 
to determine the waveguiding properties of the structure [5],[1],[6],[8], [7]. 

The structure to be analysed as shown in figure 2.1 can be divided into 
regions in which there is no change in the permittivity in the direction of 
propagation and, regions in which the permittivity changes along the direc- 
tion of propagation. 

2.1 Maxwell equation’s and Wave equation 

In the solution of any electromagnetic problem the fundamental relations 
that must be satisfied are given by Maxwell equations. 


= dB 

v * £ = -aF’ 

(2.1) 

. - dD 

(2.2) 

II 

t> 

(2.3) 

V • B = 0, 

(2.4) 


a 
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where, E is the electric field vector, H is the magnetic field vector, D is the 
electric displacement density vector, B is the magnetic flux density vector, 
J is the current density vector and p is the volume charge density. 

In addition to these relations there are three relations that concern the 
characteristics of the medium in which the fields exist. These are the consti- 
tutive relations. 


D = eE, 

(2.5) 

B — pH, 

(2.6) 

J = crE, 

(2.7) 


where, e, p, cr are the permittivity, permeability, and conductivity of the 
medium, which is assumed to be homogeneous and isotropic. For a source 
free region J = 0 and p — 0. 

Taking curl on both sides of equation (2.1) we have 


„ 0(V x H) d 2 E 

V x (V x E) = -pe- 


at 

From the vector identity V x (V x A) = V(V 
vector we can reduce equation (2.8) to 


dB 


( 2 . 8 ) 


A) — V 2 A where, A is any 


V(V • E) - V 2 F = -pe 


d 2 E 


dt 2 


(2.9) 




From equation (2.5) and (2.3) for source free region we have 


V • (eE) = 0. (2.10) 

Using vector identity V • (<M) = <?iV • A + A ■ V<£ where, A is any vector and 
4> is any scalar, equation (2.10) reduces to 

eV • E + E ■ Ve = 0. (2.11) 


From equation (2.11) it can be seen that if e is a constant (regions in 
which there is no change in permittivity in any direction), then V • E = 0, 
and equation (2.9) reduces to 

f) 2 F 

V*E - pe— = 0.. (2.12) 

This represents wave equation in uniform regions. For an assumed exp(-jut) 
time variation of the electric field, equation (2.12) reduces to 

V 2 £WV£ = 0, (2.13) 


where £ is the electric field phasor, and E(x, y , z, t) = lZe[£(x, y , z)exp(-jut)]. 

Now let e be a function of direction of propagation x, and let K(x) = ^ 
[12], [11] . Then from equation (2.11) we have 


V-E = 


E • VK(x) 
K(x) 


(2.14) 


Substituting equation (2.14) in equation (2.9) we have 


„ r E- VK(x), _2 ft Tr , ,d 2 E n 

v[ - ; ] + v‘e - = o. 


K(x) 


at 2 


(2.1S) 


For an assumed exp(—ju>t) time variation of the electric field equation (2.15) 
reduces to 

V[ g 'jr v~^] + V 2 5 + u 2 pe 0 K(x)S = 0. (2.16) 

K(x) 

Equation (2.15) represents the Helmholtz equation for non uniform regions. 



2.2 Basic equations 


To find out the waveguiding properties of the rectangular slab shown in figure 
2.1, Helmholtz equations (2.13) and (2.16) needs to be solved and the type of 
waves that may appear in every region need to be calculated. The solution 
is then constructed by choosing a suitable combination of these waves so 
as to satisfy the boundary conditions. We assume the width of the slab to 
be infinite in y-direction, and the fields are invariant with respect to the y- 
coordinate, i.e., = 0. And hence the fields can be decomposed into TE 

and TM modes, the Helmholtz equation for the field vectors then reduces 
to a scalar equation. If the layer with taps, with its width equal to the tap 
depth were not present then, the characteristic waves in the uniform regions 
appear in the form 

Si(x, z ) = £^exp[j(l3x + k z z)], (2.17) 

where, £q^ is a constant, i represents the ith uniform medium, e,- the per- 
mittivity of the corresponding medium, ft is the propagation constant in x 
direction, k. propagation constant in z direction, and k, are related by the 
equation 

,B 2 + [A4 0 ] 2 = (2.18) 

where, k 0 = y is the plane wave propagation constant in air and A the free 
space wavelength. In the absence of absorption and scattering losses, uniform 
layered structures can support surface waves that propagate along x and 
decay away from the structure in z direction. Such waves are characterised 
by real values of ft and imaginary values of k, . However, when a perturbed 
layer is superimposed on the uniform layered structure, the fields are modified 
to satisfy the boundary conditions on the taps. Under such conditions, fields 
in all. the uniform regions can be represented as a sum of plane waves, and 
are of the form, 

£i{x,z)= X £^exp[j{/3 n x + k zn z)], (2.19) 

n=—oo 

where, £$ represents the amplitudes of each partial field. The quantity /?„, 
which is the propagation constant of the nth partial field is related to the 



fundamental propagation constant /? 0 by the relation 


,3 n = Po + — , (n = 0,±1,±2,---) (2.20) 

where, d is the distance between successive taps. Each £, now satisfies the 
wave equation. Hence, kQ = ±[kf — ft 2 ]*. 

Fields in the nonuniform regions should satisfy equation (2.16). Now, 
for the case of a rectangular slab we have assumed, all the fields can be 
decomposed into TE and TM modes. For TE modes, only £ y , H x , H z are the 
only nonzero components. Now, since only y- component of the electric field 
is non-zero and also since the permittivity is a function of x the first term in 
equation (2.16) reduces to zero. The Helmholtz equation in the perturbed 
region of the waveguide can then be reduced to the scalar equation for TE 
modes as 

V% + f 2 (x)S y = 0, (2.21) 

where, f 2 (x) = u 2 pe 0 K(x) . If we assume p = p 0 (non-magnetic regions) 
then equation (2.21) reduces to 

f 2 (x) = u, Vo A'(x), = k 2 K(x) (2.22) 

where, k Q is the plane wave propagation constant in air. The electric and 
magnetic field vectors in a TE mode are related by the equations, 

£ = S y u y and H = — (V x £) (2.23) 

Now, f 2 (x) can be represented as 

f 2 (x) = k 2 p n exp(j^x) (2.24) 

n=-oo 

where p n $ are the fourier transform values for a given K(x) and represents 
the spatial frequency harmonics arising from perturbations in the cladding 
structures. Then we can assume a solution for equation (2.21) of the form 

£y= jc q n (z)exp(iPnx) (2.25) 

n=-oo 

where f3 n is given by the equation (2.20) and q n ( z ) an y function of z. 
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Substituting equations (2.24) and (2.25) into equation (2.21) we have 


Now, 


( P£ U (PS,, 

# + ^ + = 0. 


tlf n=oo 

= X ~Pl<ln{z)expUM, 

d% *-» (pq n (z) 

~TJ = X ------ - expijftx). 


(2.26) 


(2.27) 


(2.28) 


/ 2 (- r )£y = M X piexp(jlx)][ X qn(z)exp{j/3 n x)]. (2.29) 

— oo n=-oo 

Adding equations (2.27) (2.28) and (2.29) we have 

n—oo l=co 

a nq n {z)exp(j p n x) + X &/e:cp(j/3 n :r)+ 


n=oo l=oo 


A ’o X! X C'„,/ea:p(;'/? n x + i— j) = 0, 

n=-oo /=— oo 


(2.30) 


where 

ii = ^ 

and C„,/ = <?„(-)/>*• 

Now from equation (2.20) we have = fio + therefore equation 
(2.30) reduces to 

n ^° /x ,. , .2 jto . ,2x/ . 

X angn(^)exp(j/? 0 ^ + ;-7-^)+ X oze-rpO/foE + ;—£)+ 

n=-oo /=- oo 


"E E C„,,ex P (j,3 oa : + i2d!Lii) l)= 0. 


(2.31) 


n=-oo l=-oo 


Equation (2.31) is satisfied if 


k = ~[^o X CW-n-a/] 


because, the functions exp[j(/3o + ^p)x] and exp[j (/?o + 27r ^- + ^ )x] are orthog- 
onal to each other. Hence, the sum of coefficients in equation (2.31) can be 
equated to zero. Thus, for say n = 0 we have 

- -[k%(- • • +q- 3 {z)p 3 +q-i(z)p 2 +q-i{z)pi-+qo{z)po+qi(z)p-i + • • -)-pQqo(z)] 



similarly, for n = 1 we have 

~ J"T~ = ~[^o( H-2( z )P3+ < l-i( z )P2+qo{z)pi+qi(z)po+q 2 (z)p. 1 i )— j3 2 qi(z)] 

Hence, equation (2.31) reduces to a system of differential equations 

d 2 q(z) f- 

-jy- = -Pq{z), (2.32) 

where P n i = &op n -f — l3 2 8„i, P n i and 8 n i representing elements of matrix P 
and kronecker delta respectively. Since, equation (2.32) represents system of 
differential equations with constant coefficients we may assume the solution 
of the form 

q = cexp(±jaz ), (2.33) 

where c is a constant vector and a is the propagation constant in z direction. 
Substituting (2.33) back into equation (2.32) we have 

Pc - a 2 c = 0. (2.34) 

From equation (2.34) it can be seen that a 2 m is an eigenvalue of the matrix P, 
and c m represents the eigen vector for the corresponding value of a 2 m . Each 
eigen vector c m now contains elements c m „, so that every complete solution £ tJ 
in the nonuniform region now contains an infinite set of spatial frequencies, 
each of which has an amplitude c mn . Hence, for equation (2.32) we have a 
pair of solutions which can be written 


VrnK 2 ) = C m exp(-ja m z), 

(2.35) 

} ( z ) = C m exp(ja m z), 

(2.36) 


where q^\z) represent waves that travel along positive and negative z axis 
respectively. 

As cam be seen from the discussion above, the order of the matrix P 
which appears in the calculation of the fields in non uniform regions is infi- 
nite. However, in order to implement the spectral/mode-matching method 
successfully on computer for the calculation of the fundamental propagation 
constant /3 q for a rectangular slab with taps, it is necessary to truncate the 
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matrix judiciously. The matrix is of Hill’s type and hence, the truncation 
need not be unduly large for accuracies required by practical considerations. 
The truncation is found to be valid when we extend the present case to the 
theory of infinite determinants. We find such a truncation is valid provided 

[5), [16], [IS] 

(A)l*3ft>-#£-<Sl>EW. (2-37) 

K 0 i 

which must hold for |n| > N where N is a finite positive integer and the prime 
in the summation indicates that the i = n term is excluded, ko is the plane 
wave propagation constant in free space, and p n s indicates the value of the 
Fourier transform for a given profile. 

2.3 Field solution for a single non uniform 
layer 

In this section we discuss the problem of calculating the fields supported by 
a single non uniform layer surrounded by uniform regions as shown in figure 
2.2. This solution can then be extended to the actual case of a rectangular 
slab waveguide with taps in cladding. This is discussed in the next section. 



14 


In the figure 2.2 region 1 represents either core or cladding depending on 
the tap depth, region 2 represents air, and region 3 is the region with taps 
i.e., the region in which refractive index varies in the direction of propaga- 
tion. A plane wave of unit amplitude is incident from region 1 and it sets 
up fields in all the regions. The fields that are setup in various regions can 
be obtained from the discussion carried out in previous sections this chap- 
ter. Using equations (2.25) (2.35) and (2.36) the electric and magnetic field 
components within the nonuniform layer can be written as 

m=o o n=oo 

£ 3 - J2 [g^ expti <*mz) + g^expi-j a rnz)] V mn exp(j /3 n x) , (2.38) 

m=~oo rt=— 00 

m= 00 n= 00 

7*3= X2 [g^expijQmz) -g { ~ ] exp{-ja m z)] 52 I mn txp{jfi n x), (2.39) 

m=—oo n=—oo 

where, c mn = V mn for TE modes and c mn = i mn for TM modes. If one set is 
known, other can be calculated using Maxwell equations, g W and gj~^ are 
unknown constants. Using equation (2.19) the fields in region 1, which is a 
uniform region can then be given by 

rc=oo 

£1 = exp{j/3 0 x + jkil ) z)+ 52 s n exp({jp n x - jk£}z), (2.40) 

n=~ 00 


Hi = Y^ l) exp(jp 0 x + j&UM - 51 Y^SnCxpdjpnX - jk^z), (2.41) 


where, s n are the amplitudes of the scattered waves. Similarly, fields in region 
2 can be given by 


71=00 

£2 = 5Z *»» eX 

n=— 00 


p{jPnX + jk { „ 2 Jz)i 


(2.42) 


= 52 tnY^expijpnX + jk$z), 

n=— 00 

where, i s the characteristic modal admittance given by 



jfe(0 

— £2- for TE modes and , 
ojpo 

= t-pr- for TM modes 
u(') 

fczn 


(2.43) 

(2.44) 

(2.45) 



and t n are the amplitudes of radiated waves. Now, applying the boundary 
conditions at z = 0 for the continuity of electric and magnetic fields and 
using equations (2.38) (2.39) (2.40) and (2.41) we have 


m=oo 


■5o» + s»= E K.nbi. +) + £■>]. 

m= — co 

(2.46) 

771=00 

m=-oo 

(2.47) 


similarly, applying boundary conditions to at z = t g to equations (2.38) (2.39) 
(2.42) and (2.43) we have 

771 = 00 

£ VmniexpijcimtgjgW + exp(-ja m t g )g ^ J = exp(jkj$t a )t n , (2.48) 

m=— oo 

m= oo 

£ I mn [exp(ja m t g )g^ ) - exp{-ja m t g )g^ ) ] = exp(jk<$t g )Yi 3) t n . (2.49) 

771 = — OO 

The equations (2.46) to (2.49) can be written in matrix form, the order 
of the matrices is equal to the order of matrix P. The matrix P is truncated 
using equation (2.37). If N is the order of the truncated matrix then, equation 
(2.46) can be written in matrix form as 


e + s = V]<7 (+) + <7 ( ’] 


(2.50) 


where 


e = 


/0\ 

1 

\0j 


t S == 


f S=N \ 
2 


So 


\ / 


( ff-jl \ 


J-' = 


9o 


\ 4 ’ / 


,5 (+) 


and elements of matrix V = {c mn } for TE modes. 
Similarly equation (2.47) can be written as 

f 'i(e-S) =/[j (+l 


/attJ \ 
2 


9 0 


v 4 +) j 


(2.51) 


where Y\ = diag{l^^}, and I = {cm,,} for TM modes. 



Equation (2.48) can be written as 


V[C’ g m + S W ] = Eh i (2.52) 

where = diag {exp(ja m t g )}, ~E' rp > = diag {exp(-ja m t g )} and 
~E kz = diag {exp(jk^}tg)}. similarly equation (2.49) can be written as 

9 W - Sip’ j H l = Yi%. t (2.53) 

where Y% = diag{Fj 2 ^}. 

Now, premultiplying equation (2.52) with F 2 and subtracting from equa- 
tion (2.53) we have 

(I - Y2V02 j< +) - (I + Y WjSipjM = 0 (2.54) 

From the above equation we have 

s<-> = 1 » (+i 

where = (/ + V' 2 F) _1 (/ - f' 2 V). Now, since is a diagonal matrix 

=( * 4 ") 

its inverse is equal to E xp , hence the above equation reduces to 

a<-> = lll’ipi 1V +I (2.55) 

Pre multiplying equation (2.50) with F i and adding to equation (2.51) we 
have 

2F x e = (/ + YiV)g {+) -(I- Yifyg^ 
substituting for from equation (2.55) in the above equation we have 

2Y t e = ((/ +Y l V)-(I- K, V)(f * ) R,llp ) )]a< + > 

Pre multiplying the above equation with (I + F x V’) -1 we have 

S 9 g + = foe, (2.56) 

where 

S, = I - (2-57) 

T 0 = 2(/ + fif , )-‘? l , (2.58) 
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R 0 = (I+Y 1 V)-\I-f l V). (2.59) 

I is the identity matrix, To and Rq are transmission and reflection matrices 
looking into region 1 at z = 0. 

In the absence of any incident wave from an exterior region, we have null 
vector instead of e in equation (2.56) which is satisfied only if the determinant 
of S g vanishes. Therefore 

zs — ~ — 4 - ) — .■■ ( 4 - ^ 

det(Sg) = detfl - R Q E xp R g E xp } = 0. (2.60) 

Equation (2.60) represents dispersion relation, solving which yields the un- 
known value of the propagation constant /3 0 . An iterative approach has 
been followed to solve for f3 0 which is described in the next chapter. Sim- 
ilarly by solving for g ^ and g^ the fields in the tapped regions can be 
calculated, and by solving for s and i the scattered and radiated fields can 
also be calculated. 

2.4 Field solution for rectangular slab with 
rectangular taps 

The aim of our thesis is to find out the wave guiding properties of the a rect- 
angular slab with taps in cladding as shown in figure 2.1. The problem can 
be easily solved by making a simple extension to the discussion carried out 
in the previous section. The rectangular slab with taps in the cladding can 
be treated as superposition of additional uniform layers below the structure 
shown in figure 2.2. The problem is further simplified by recognising the fact 
that superposition of uniform layers below the non-uniform region can be 
treated as change in admittance matrix at the plane z = 0. And also from 
the fact that a single non uniform layer can be represented by equivalent 
transverse network shown in figure 2.3. 

In this network, each semi-infinite transmission line (shown by thick line 
in figure 2.3) represents one of the modes in each region. All the transmission 
lines are connected to the non uniform region described by the box which 
represents the matrix S g . Now the additional uniform layers can be added 



below by simply interpolating set of transmission lines of appropriate lengths, 
as shown in figure 2.4. 

In this case at the plane z = —tf we define reflection coefficient 


yd 1 ) _ yd3) 

there by the input admittance at z = 0 is given by 

yin _ y(l) l ~ Pnexp[2jk^tf) 

1 + p n ex p(2jkintf) 


(2.61) 


(2.62) 


Hence the reflection coefficient matrix R 0 as described by equation (2.59) 
gets modified to 

Ro = (I+ Yin V)- 1 (I - Yin V). (2.63) 

There by the matrix S g which is dependent on R 0 also gets modified and 
hence the propagation constant. Hence to find out the propagation charac- 
teristics of a rectangular slab with rectangular taps it is necessary to make 
the determinant of the modified matrix S g zero. 


2.5 Field solution for rectangular slab with 
curved profile taps 

In the discussion till the previous section we have concentrated on rectangular 
taps, however using the spectral/mode-matching method analysis can be 
done for taps with curved boundaries as well. The approach is to judiciously 
approximate the curved profiles. Analysis has been carried out for taps 
having triangular and gausssian shapes as well, as shown in figure 2.5. 

The approximation can be done by partitioning the layer with taps into 
fine layers and approximating each of these by rectangular profile. We thus 
obtain a configuration as shown in figure 2.6. The analysis of such multiple 
layered structure can be carried out by straight forward extension of the 
discussion carried out in the previous sections. It can be seen that the matrix 
Rg is the reflection coefficient matrix looking into region 2 at z = t g in figure 
2.2. When a curved profile is approximated by rectangular layers, the matrix 
R g gets modified and hence the propagation constant. 



A much simpler, but less accurate [5], way for analysing a curved profile 
is to average the permittivity for every value of z inside the region with 
taps. Thus, if the curved boundary that seperates two media that have 
permittivities ej and e 2 is described by a function 

z = h(x ), for 0 < z < t g 

the averaged permittivity then becomes 

e(z) = e 2 + (e 2 - 

*9 

The problem then reduces to that of a layer with uniform thickness but, 
varying e(z) and can be solved easily. Since, now the curved profile is ap- 
proximated just as a change in permittivity in x-direction and hence the 
problem is similar to that of a rectangular profile, but with a continuous 
change in the permittivity. 

From the discussion that has been carried out in this chapter it is clear 
that the method can be applied to practically any dielectric waveguide. A 
simulation package based on the spectral/mode-matching method that has 
been described in this chapter, has been developed. The package is described 
in the next chapter. 







Figure 2.5.2 : Gaussian taps 
Figure 2.5: Rectangular slab with curved profile taps 
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Chapter 3 

Simulation package for 
spectral / mode-matching 
method 


Based on the spectral/mode-matching technique which was described in the 
last chapter, we have developed software which gives the propagation char- 
acteristics of rectangular slabs with taps in cladding. In this chapter we 
describe the software that has been developed. We also discuss the validity 
of the software. To test the validity of the software, propagation characteris- 
tics of a rectangular dielectric waveguide were calculated theoretically. The 
propagation characteristics of the slab were then calculated by the software 
by making the tap depth equal to zero. The error was found to be less than 
1% of the theoretical value at all wavelengths.. 

3.1 Various issues in software implementa- 
tion 

The application of the spectral/mode-matching technique described in pre- 
vious involves the solution of the transcendental equation, which involves 
setting determinant of S g given by equation (2.60) to zero. All the parame- 
ters entering S g axe known except for /?o which is regarded as the unknown 
variable for any given frequency. In general however, the quantities a m and 
c mn (in the region with taps) which are the eigen values and eigen vectors of 
matrix P axe not known explicitly, so that determination of their values is 
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paxt of the computation process for finding j3 0 . 

3.1.1 Truncation of the matrix of infinite order 

Truncation of the matrix of infinite order for the calculation of the eigen 
values of matrix P forms one of the most important parts of the computation 
process. The matrix is of Hill’s type and hence the matrix need not be of 
unduly large order. The validity of truncation is checked by making use of the 
equation (2.37) [16], [18], [5]. To start with, the eigen values are calculated 
for an assumed order of the matrix, say five. The condition for validity is 
checked, if the condition is not satisfied the order of matrix is increased by 
one. In general however the order of the matrices to be dealt with was found 
to be sufficiently small, typical values being 9 to 15. However, the increase 
in the order of the matrix increases the accuracy of the result, at the same 
time results in increase in computational time. 

3.1.2 Basic structure of the subroutines to be imple- 
mented 

A first step in the programming of a computer routine for solving the tran- 
scendental equation described by S g is to provide routines for calculation of 
the Fourier transform given the tap shape, width, distance between succes- 
sive taps number of taps, the refractive indices of cladding and tap filling if 
any. Since for the tap shapes we were dealing with the Fourier transform can 
be derived analytically, we used the derived results in calculating the Fourier 
transform. Once the Fourier transform has been calculated the next step is 
to find out matrix P for an assumed value of /3 0 , then the eigen values a£, 

, the eigen vectors c m for an assumed order. The validity of truncation is 
checked. After the truncation is found to be valid, the matrix S g is found, 
the determinant is calculated. If the determinant is found to be equal to 
zero, the value of /?o is that of the actual propagation constant. However, 
if the determinant is found to be not equal to zero, then the value of flo is 
varied till the determinant of S g is found to be equal to zero. 

All the routines have been written in ”C++” programming language, 



except for the routine which calculates the eigen values and eigen vectors of 
a complex matrix of any order. For calculating the eigen values and eigen 
vectors of a complex matrix NAG routine has been used, the routine uses 
similarity transformations to reduce the given matrix into Hessenberg form 
and then calculates the eigen values and eigen vectors. The flow chart that 
has been followed is given in figures 3.2 and 3.3. 

3.1.3 The main loop 

After reading all the input parameters, the program calculates the Fourier 
transform values for a given set of tap parameters. Since, for the tap shapes 
we were dealing with the Fourier transform can be derived analytically, we 
used the derived results in calculating the Fourier transform. The program 
then enters the main loop (shown in figure 3.2). The range of wavelengths 
over which we have calculated the propagation characteristics is 0 .8fim to 
1.6//m(A max ). This range was chosen since, it is the range of interest for 
optical communication. The program calculates the propagation constant 
for each wavelength in this range starting from 0.8 fim with an increment of 

0. 01/im(A, ucr ). The basic idea behind the program is to calculate a suitable 
value of d 0 which would result in making the determinant of matrix S g zero. 
In order to do this we start with an initial guess for do, which was taken to 
be equal to 1.01 where, n 2 is the refractive index of cladding, A being the 
wavelength of operation. This value was chosen since it is known that for 
a rectangular slab the value of propagation constant lies between pp and 
pp. Using the initial guess and the fourier transform values the elements of 
matrix P are found for an assumed order. The eigen values and eigen vectors 
of the matrix are then calculated using NAG routine. Since the matrix P 
is truncated, condition for the validity of the truncation is checked. If the 
truncation is found to be invalid then, the order of the matrix is increased by 

1. If the truncation is found to be valid then the program proceeds further 
and the elements of matrix S g are calculated using the equations derived in 
the previous chapter. The determinant of matrix S g is then calculated. If 
the determinant is found to be greater than 10 -6 , the value of determinant 



is stored in some temporary memory location, the value of /3 0 is increased by 
5%, the determinant of matrix S g is again calculated following all the steps 
described above. If the value of determinant turns out to be less than 10 -6 
then the value of 0o gives the value of the propagation constant. On the other 
hand if the value determinant is again greater than 10 -6 , the present value 
of the determinant is compared with the previous stored value, if the present 
value of determinant is less than the previous stored value the value of /3 0 
is further increased otherwise, the value of /? 0 is decreased. The method is 
thus iterative and finally gives the value of 0 O , which makes the determinant 
of matrix S g zero (less than 10~ 6 ). The value of A and its corresponding 
value of fio are written in a output file. The value of A is now incremented by 
Xincr the value of /?o is calculated by employing the method described above. 
This is repeated until the value of A reaches 1.6 /j,m(\ max ), after which the 
program comes out of the main loop. The program thus calculates the value 
of the propagation constant iteratively. 

3.2 The simulation program 

In this section we describe the simulation program that has been written in 
"C++” for the implementation of the spectral/mode-matching method . 
Input parameters : 

1. The software developed can handle three tap shapes rectangular, tri- 
angular, and gaussian. 

2. In case of rectangular and triangular (base width of triangle) taps the 
width of tap (dl) is the input, however in case of gaussian taps, the 
width of the tap has been taken to be 10 times variance of the gaussian 
profile (shown in figure 2.5). 

3. Number of taps (N). 

4. Distance between successive taps (d). 

5. Refractive indices of core (nj), cladding ( 712 ) and refractive index of 
tapfilling material if any. 
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6 . Total width of the slab (D). 

7. Width of the core (t/). 

S. Depth of taps (t 5 ). 

9. Name of the output files. 

Output parameters : 

The program stores the output data in two files one of which contains 
wavelength vs propagation constant data, and the other contains the error 
in making the determinant of matrix S g zero. 

3.3 Validity of the simulation program 

To check the validity of the program that has been written, we decided to find 
out the propagation constant at various frequencies by making the tap depth 
equal to zero in our program , and then compare the values so obtained, with 
the values of propagation constant obtained by theoretical calculation for a 
rectangular slab. 

3.3.1 Theoretical solution 

The electric field of the fundamental TE mode[14] for a planar waveguide is 
given by 

E y (z) = >lcos(/c 2 :)for|x| < — (3.1) 

and 

E v {z) = Se^for|x| > ^ (3.2) 

where d is the core width. A and B are found by equating Acos(kx) and 
at x= | . k and 7 are given by 


II 

(3.3) 

2 v 2 0 

7 = 3V T~ e 

(3.4) 



where V is given by 


V = ~ nl (3.5) 

where n l5 n 2 are refractive indices of core and cladding respectively and A 
the wavelength, f is obtained as the solution of the transcendental equation 

jV * 

^tan(0 = y— (3.6) 

This transcendental equation was solved by the bisection method. The lower 
and upper limits of the propagation constant being and 2 h*I. The 
propagation constant was obtained for A varying from 0.8 [xm to 1.6^m, for ni 
= 1.4616, n 2 = 1.4571 , the width of core = 3.8 fim. For the same parameters 
the tap depth was made equal to zero in the simulation program and the 
propagation constant was thus obtained. The figure 3.4 shows the plots of 
the A vs (3 0 (propagation constant). The plot of A vs percentage error is 
shown in figure 3.5 . It can be seen that the error is less than 1% at all the 
wavelengths. 

3.4 Conclusions 

In this chapter we have presented the methodology that has been adopted to 
calculate the effect of taps on the propagation characteristics of the optical 
fiber. We also showed that the method the sufficiently reliable, and that 
the method can be applied to any dielectric waveguide, since all the input 
parameters entering the program can be varied. 
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Figure 3.3: Plots showing 











Chapter 4 

Simulation results 


In this chapter we discuss the results that have been obtained from the 
simulation experiments carried out, based on the spectral/mode-matching 
approach described in the previous chapters. The results have been obtained 
for the case of a planar waveguide whose total width and the width of the 
core were taken to be equal to the diameters of cladding and core respectively 
of a standard optical fiber. The various input parameters which have been 
kept constant through out the simulation are 

1. Refractive indices of core and the cladding are taken as 1.4616 and 
1.4571 respectively. 

2. Total width (in /.-direction) of the planar waveguide has been kept 
constant at 125 /zm . 

3. Width of the core has been taken to be 3.8 pm . 

The data is taken from a single mode fiber commercially available [19]. For 
these parameters the waveguide remains in the single mode over the A range 
from 0.8/zm to 1.6/zm. The simulation program however, is sufficiently gen- 
eral to handle any values for the above parameters. 

The parameters that have been varied are 

1. Depth of the taps ( t g ), for which we took typical values between 0 to 
60.6 /zm. For t g — 60.6/zm the tap touches the core cladding boundary. 

2. Shape of tap, which can be made gaussian, rectangular, and triangular. 
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Tap parameter j 

Typical values 

Depth j 

0.0 m to €0.6 pm 

Shape j 

Rectangular, Triangular, Gaussian 

Number 

10 to 200 

Width 

20^m to 200 pm 

Distance 

100 pm to 1mm 


Figure 4.1: Table showing typical input parameters 


3. Number of taps (N), for which typical values were taken between 10 to 

200 . 

4. The tap width (dl) was varied from 20 /im to 200pm. 

5. The distance (d) between taps was varied from 100pm to 1mm. 

4.1 Variation of tap depth 

We executed our program several times for different tap parameters. In each 
case only one parameter was varied keeping all the other parameters constant. 
Figure 4.1 shows the table of typical values for all the parameters. For each 
variation we calculated wavelength (A) vs propagation constant (/Sq). The 
graphs have been plotted with wavelength (A) on the x-axis and propagation 
constant (,d o) on the y-axis. Range for wavelength was selected to be 0.8pm 
to 1.6 pm, since it is the range for optical communication systems. 

In this section we present the results that have been obtained by varying 
the tap depth keeping all other parameters constant. Number of taps was 
taken to, be equal to 10. Tap width was taken to be equal to 20pm, in case 
of gaussian taps the width of tap was taken to be equal to be 10 times of the 
variance. The Distance between successive taps was kept constant at 1mm. 
The results of the tap depth variation are shown in figures 4.2 and 4.3. The 
figures show the variation of tap depth for the case of triangular taps. From 



the figures 4.2 and 4.3, it can be seen that the propagation constant increases 
as the tap depth is increased. This can be explained, since an increase in tap 
depth results in increased radiation from the taps, and hence an increase in 
tap depth affects the signal propagating inside the fiber However, the effect 
is more profound at low wavelengths. The propagation constant does not 
show variations above 5% with respect to the no tap case in the range 1.3 
to 1.6/mi, which is range* of interest for optical communication. Hence, it 
can be concluded that the for small number of taps (10) and tap depths 
below 40 ftm, the the* increase iu propagation constant with the increase in 
tap depth is below 5%. 

4.2 Variation of tap shape 

In this section we present the results of variation in the shape of the tap 
keeping all the parameters except the tap shape constant. The results are 
shown in figure 4.4. From tilt* figure 4.4 it can be concluded that the variation 
in tap shape does not affect the propagation characteristics significantly. It 
can also bo seen from the figure 4.4, that the triangular and gaussian taps 
have similar effect on the propagation characteristics, and hence it can be 
concluded that the gaussian tap can be approximated to a triangular tap 
for analysis. It can also be supported by the fact that the spatial frequency 
spectrum due to the triangular and gaussian taps is similar as shown in figure 
4.13. 

4.3 Variation of number of taps 

In this section we present the results that have been obtained by varying the 
number of taps keeping all the other parameters constant. The tap width 
was kept constant at 20 /im, the distance between taps was kept constant 
at 1mm, the tap depth was kept constant at 30 gm . From figure 4.5 it 
can be seen that the propagation constant increases as the number of taps 
is increased, however again the effect is more profound at low wavelengths. 
The propagation constant increases for larger wavelengths 1.3 nm to 1.6 fim 



as the number of taps is increased beyond 100. This can be explained from 
the fact that as the number of taps is increased the amount of optical energy 
that is radiated out of the taps also increases, and hence the number of 
taps have a profound effect on the propagation characteristics of the fiber. 
Figures 4.6 and 4.7 show the results of variation of number of taps in case 
of rectangular taps with tap depth equal to 50/zm, distance between taps 
equal to 1mm and width of taps equal to 20/zm. As can be seen from the 
results the increase in the tap depth from 30pm to 50/zm has more effect on 
the propagation constant even the range 1.3pm to 1.6pm. This is supported 
by the fact, that increase in number of taps results higher amplitudes of the 
spatial frequencies and the increase in tap depth results in more radiation. 
The spatial frequency spectrum showing the increase in amplitude with the 
increase in number of taps is shown in figure 4.8. 

4.4 Variation of tap width 

In this section we present the results that have been obtained by varying the 
tap width keeping all other parameters constant. The number of taps was 
kept constant at 10, the distance between taps was kept constant at 1mm and 
the tap depth was 30pm . From the figures 4.9, 4.10 and 4.11 it can be seen 
that that the propagation characteristics are not affected by the variation 
in the tap width in case of rectangular taps however, in case of triangular 
and gaussian taps the change in tap width does affect the propagation char- 
acteristics. This can be explained from the fact that in case of triangular 
taps a change in tap width keeping the tap depth constant effectively means 
the change in the angle of the facets of the tap and hence the amount of 
radiated energy also varies. In case of rectangular taps however, the change 
in tap width does not effect the angle at which the propagating mode strikes 
he facets of the taps. Since we were dealing with tap widths much larger 
ian the wavelength, the change in the width of a rectangular tap does not 
ow a large effect on the propagation characteristics. In case of gaussian 
a change in tap width means a change in the shape itself, and hence the 
^gation characteristics are effected by the change in tap width in case 
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of gaussian taps. However, even in cast* of triangular and gaussian taps the 
wavelength range 1.3 fun to l.6/wi remains unaffected. 

4.5 Variation of distance between taps 

In this section we present results that have been obtained by varying the 
distance between taps keeping all other parameters constant. The n um ber of 
taps was kept constant at 10, the tap width was kept constant at 50 gra and 
the tap depth was kept, constant at 50/im . From the figure 4.12 it can be 
seen that the propagation characteristics are not effected by the variation in 
distance between taps. This can again be explained from the fact that the 
all the distances involved are large compared to the wavelengths. 


4.6 Conclusions 


It can be concluded from the results of the simulation experiments, that the 
tap depth and the number of taps have a profound effect on the propagation 
characteristics of the optical fiber, where as the tap width and the distance 
between taps do not effect tin* propagation characteristics of the optical fiber. 
This fact can be put to use in the making of fiber optic delayline since, in 
a fiber optic deiayline the distances involved are in the range of hundreds 
of microns and the delayline should not effect the signal propagating inside 
the fiber. On the other hand the tap depth and the number of taps, which 
have a profound effect on the propagation characteristics can be used in the 
making of dispersion less and dispersion flattened fibers. 

From the results so obtained, the waveguide dispersion at any desired 
frequency can be calculated by making use of the formula [11]. 


o 

" " 1 Ac ' iV 2 

where, A is the wavelength of operation, c is the velocity of light in free space, 
rii and n-i are refractive indices of core and cladding respectively, V is the 
normalised frequency given by 


V = - nj 
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N = 1 0, Tap width = 20 pm, Distance between taps = 1 mm 



A. (p m) 

N = 10, Tap width = 20 pm, Distance between taps = 1 mm 

Figure 4.3: Tap depth variation in triangular taps 
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N = 10, Tap width = 20 pm, Tap depth =30 pm, Distance between taps = 1mm 

Figure 4.4: Comparison of tap depth variation and tap shape 







X (n m) 


Tap depth = 30 n m, Tap width =20 p. m, Distance between taps = 1 mm 

Figure 4.5: Variation of number of taps in rectangular taps 






X(p m) 



Tap depth = 50 p m, Tap width = 20 |i m, Distance between taps = 1mm 


Figure 4.6: Variation of number of taps in rectangular taps 
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Tap depth = 50 pm, Tap width = 20 pm, Distance between taps = 1mm 



Tap depth = 50 p m, Tap width = 20 pm. Distance between taps = 1mm 

Figure 4.7: Variation of number of taps in rectangular taps 





Amplitude 
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Spatial frequency 3 *“ 

Distance between taps = 1mm , Tap width = 20 pm 


Figure 4.8: Spatial frequency spectrum for variation in number of taps 





N = 10, Tap depth = 10 p m, Distance between taps = 1mm 



N = 10, Tap depth = 10 pm, Distance between taps = 1mm 

Figure 4.9: Variation of width of taps in rectangular taps 











N = 10, Tap depth = 30 p. m, Distance between taps = 1mm 

Figure 4.11: Variation of width of taps in gaussian taps 









Amplitude 



N = 5, Tap width = 20 4 m, Distance between taps = 1mm 
Figure 4.13: Spatial frequency spectrum of Triangular and Gaussian taps 




Chapter 5 
Conclusions 


In this thesis we carried out a numerical simulation of the effect of a series 
of radiative LACE made taps on the propagation characteristics of optical 
fibers. We use a spectral/mode-matching method for our analysis. The 
various parameters that have been looked into are the effect of tap shape, 
tap depth, tap width and the number of taps on the propagation constant 
/?o. The following conclusions were obtained from our observations of the 
results that we obtained. 

• The effect of taps is to increase the propagation constant. 

• As the number of taps is increased the propagation constant also in- 
creases. 

• The propagation constant increases with the increase in the tap depth. 

• The propagation constant increases with the increase in the width of 
the tap, however not much increase was observed when compared with 
the other parameters. 

• The use of the package so developed to investigate the propagation 
characteristics of any dielectric waveguide is also evident since, all the 
parameters entering as inputs to the program can be varied. 

• The results so obtained show a relatively strong dependence of the taps, 
on the propagation characteristics of the optical fibers, when compared 
with the other parameters. 
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5.1 Suggestions for future work 

The following are the suggestions for the further work that can be done, to 
understand more about the behaviour of tapped optical fibers. 

• The investigations have been carried out on the rectangular slab waveg- 
uide in this thesis, however the method can be extended to the actual 
optical fiber. 

• Further investigations can be carried out on multi-mode fibers and 
rectangular slabs as well. 

• In this thesis we have concentrated on the general solution of the 
Helmholtz equation, further work can be carried out to obtain the 
particular solution, so that radiative properties of the taps can also be 
studied. 

• The effect of the material dispersion on the propagation characteristics 
can be included. 

• A similar study can be carried out by filling the taps with various 
dielectric materials. 

• The seemingly strong dependance of the various tap parameters on 
the propagation characteristics of the optical fiber, can be used in the 
design of optical fibers with the desired value of propagation constant 
at a given frequency. This property can then be utilised in design of 
dispersion flattened and dispersion shifted fibers. 

To conclude it can be said that the laser etched fiber optic taps do have 
a significant effect on the propagation characteristics of the optical fibers. It 
can be expected that they play a important role in future in optical fiber 
communication systems, and this work is only a small step in that direction. 



Appendix A 
Fourier transform 


The spatial fourier transform of a function f(x) is given by 

/ CO 

f{x)^~ jkx Ux (A.l) 

-oo 

where k is the spatial frequency. 


dl d 

-<= =». 
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x=0 

Figure A.l: Profile for Rectangular taps 


The fourier transform for the rectangular tapped structure can be derived 
by viewing the taps as a change in permittivity of the material, ei is the 
permittivity of tap and e 2 the permittivty of the cladding. The taps are 
assumed to be made from x=0 to x=Nd, where N is the number of taps and 
d is the distance between taps. If the width of tap is assumed to be equal to 
dl. The fourier transform for a single tap is then given by 

F(k) = f ^ e^-^dx + I* e 2 e ( ~ jkx) dx (A.2) 

Jo J dl 

by integrating, we have 

F(k ) = -rr[ei(l - e^ kd1 ^) + e 2 {e^ jkd ^ - e^' w >)] (A.3) 

jk 
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Amplitude 



Figure A. 2: Fourier transform 


The fourier transform for the entire structure is then given 

G{k) = F(k)[l + e(' jd) + e ( - jM) + • • • + e^ N ~ 1)d) } 

the above equation then reduces to 

1 _ e (-jNd) 

m = F(k)[ ■ 


(A.4) 


(A.5) 
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